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We study the dynamical stability of "tachyacoustic" cosmological models, in which primordial 
perturbations are generated by a shrinking sound horizon during a period of decelerating expansion. 
Such models represent a potential alternative to inflationary cosmology, but the phase-space be- 
havior of tachyacoustic solutions has not previously been investigated. We numerically evaluate the 
dynamics of two non-canonical Lagrangians, a cuscuton-like Lagrangian and a Dirac-Born-Infeld La- 
grangian, which generate a scale-invariant spectrum of perturbations. We show that the power-law 
background solutions in both cases are dynamical attractors. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The physics of the very early universe is a rich arena 
for theory, and there are an abundance of potential mod- 
els that solve, at least partially, the well-known problems 
of the standard cosmological paradigm. Simple single- 
field inflation 1] is surely the most successful model, but 
alternatives to it have been proposed: pre-big bang cos- 
mology [2], ekpyrotic and cyclic models [Sj, nonsingular 
quantum cosmological models [1], non-canonical models 
[5], among others. It was recently shown that in an ex- 
panding cosmology with standard General Relativity, it 
is possible to generate a spectrum of super-Hubble cur- 
vature perturbations consistent with data in only three 
ways [6]: 

• Accelerating expansion (i.e. inflation), 

• A speed of sound faster than the speed of light, 

• Super-Planckian energy density. 

In this paper we focus on the second of these possibili- 
ties, a superluminal sound speed. In 7J we investigated a 
method of solving the cosmological horizon problem and 
seeding scale-invariant primordial perturbations in a cos- 
mology with decelerating expansion and a corresponding 
growing comoving Hubble horizon. If one has a decay- 
ing, superluminal sound speed, curvature perturbations 
can be generated outside the Hubble horizon without in- 
flation. We proposed the term tachyacoustic for such 
cosmologies, which are closely related to varying speed 
of light theories. 

Using a generalization of the inflationary flow formal- 
ism [8 introduced by Bean, et al. [9-llJ for noncanon- 
ical models, we derive two Lagrangians with solutions 
which exhibit a shrinking comoving sound horizon and 
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decelerating expansion. The first is a "cuscuton" La- 
grangian |12| . which is linear in the scalar field kinetic 
term, instead of quadratic as in the case of a canonical 
Lagrangian. The second is a Dirac-Born-Infeld (DBI) 
Lagrangian, of a form similar to that used to implement 
DBI inflation in string theory [T31 [H] • 



The cuscuton Lagrangian is particularly interesting 
from an observational perspective. In a recent paper |15j 
it was shown that the cuscuton-like model also leaves 
a potentially observable non-Gaussian signature in the 
CMB anisotropy field, /nl ~ which is slightly dif- 

ferent from the signatures of other superluminal models 
discussed in the literature [THl E] . This difference is due 
to the contribution of the term C,^ in the cubic action 
[ISl [12] to the full non-Gaussian amplitude, giving rise 
to an extra term with a linear dependence on the flow pa- 
rameters e and s. This dependence is not present either in 
DBI models (since the coefficient of C,^ in the third-order 
action is identically null in this case), or in disformal bi- 
metric model |16| . since by projecting its scalar- field ac- 
tion in the Einstein frame one obtains a DBI-like action. 
Therefore, the cuscuton tachyacoustic model "inherits" 
an extra contribution from the C,^ term, providing a dif- 
ferent value for /nl compared to the other superluminal 
models. Although small, this difference might be observ- 
able in the near future, which is of great importance for 
falsifying superluminal noncanonical models. 



However, a detailed study of the stability of the dy- 
namical system for these tachyacoustic models is lacking, 
and we fill this gap in the present paper. We investigate 
the dynamics of both the cuscuton and DBI Lagrangian, 
and show that solutions in both cases are dynamical at- 
tractors. This paper is organized as follows. In Section 
we review the basics of tachyacoustic cosmology. In Sec. 
|III|we discuss the attractor behavior of the cuscuton-like 
and the DBI model (IIIIBl). In Sec 
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II. TACHYACOUSTIC COSMOLOGY 

Tachyacoustic cosmology is a particular solution of a 
wider class of k-essence models (see Bean et al. [S]), 
which we briefly review below. Consider a general 
Lagrangian of the form C = C [X, (j)] , where 2X = 
9^^ dii4>du'i' is the canonical kinetic term (X > accord- 
ing to our choice of the metric signature). The energy 
density p and pressure p are given by 



p = 2X£x - 
The speed of sound is given by 



PX \ L-X 



(1) 

(2) 



(3) 



and the corresponding equation of motion for the field 
is 



C 



X 



0, 



(4) 



where the subscript "X" indicates a derivative with re- 
spect to the kinetic term, and the subscript indicates a 
derivative with respect to the field. The Hubble param- 
eter H is determined by the Friedmann equation, 



1 



1 



3M|,' 



{2XCx - C) , 



and the continuity equation is 



p ^ 2HH = -3H {p + p) = -6HXC 



X, 



(5) 



(6) 



where we use the the reduced Planck mass Mp = 
I/VSttG. For monotonic field evolution, the field value 
4> can be used as a "clock" , and all other quantities 
expressed as functions of <j)^ for example X ~ X{(j)), 
C = £[X {(f>) ,(/)], and so on. We consider the homo- 
geneous case, so that = \f2X. Next, using 



dt 



(7) 



we can re-write the Friedmann and continuity equations 
as the Hamilton Jacobi equations. 



/2X 



L-X 



C 



C. 



X 



(8) 
(9) 



where a prime denotes a derivative with respect to the 
field (p. The number of e-folds N is defined as the loga- 
rithm of the scale factor, 



a cx e 



-N 



(10) 



and can be re-written in terms of dq 



dN = -Hdt = 



'2X 



2Ml 



by: 

H'i<f>) 



dcj). (11) 



As in the case of canonical inflation |8], we can in- 
troduce a hierarchy of flow parameters, the flrst three 
defined by 



1 dH 



2Ml 



H'{4>) 



HdN Lx{(f)\H{4>) 

1 dcs _ 2Ml H' {(j)) c's {(j)) 
'^dN~~Cx (0) i?(</>) cs ((/.)' 
1 dCx _ 2M2 H' (0) C'x (0) 



Cx dN £^(0) H{<l>) £x(0) 



(12) 
(13) 
(14) 



We construct an exact solution such that the parameters 
e, s, and s are all identically constant, the background 
evolution is a power law, 

a cx e"^ cx t^/\ (15) 

and the Hubble parameter and sound speed evolve as 



cx e 



eN 



H 



cs cx e 



X cx e 



sN 
sN 



(16) 



In terms of the field 0, this solution to the flow equations 
corresponds to 



Cx = 



8e fMp 



H (0) = Ho 



cs (0) 



00 



-2c/s 



2s/ s 



(17) 



(18) 



(19) 



where 0o is a fiducial field value defined such that 
Cs (</'o) = li and the field evolves as 



-sN/2 



(20) 



From Eqs. ( |14|17[ ), we see that 0o is not an independent 
constant, but is given by 



Mf 



Be 

22- 



(21) 



From the Hamilton- Jacobi Equation (|8| , the field velocity 
is then 



2 \ (PQ 



l-2e/s 



(22) 
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We then construct a Lagrangian for which Eqs. 
( 17|18|19 ) are an exact solution to the associated equa- 
tions of motion. The important point is that given a 
solution to the flow equations, the associated Lagrangian 
is fully determined, up to the specification of a choice of 
gauge, i.e. the relationship between s and s [9 . For a 
particular choice of gauge, we construct the Lagrangian 
as follows 7\: From Eqs. (17) and (19 1, we see that the 



speed of sound eg can be written in terms of Cx 



4 



1 + 2X'- 



-I -1 



XX 



C 



X 



where we have used Eq. ([s]), and defined 

_ ( 



C 



(23) 



(24) 



The result is a differential equation for the Lagrangian 
£(X,</.): 



2XCxx +Cx-CCl=0 



where we have defined 



1 + ^. 

s 



(25) 



(26) 



Therefore, by specifying a relationship between the pa- 
rameters s and s, we can construct a Lagrangian as the 
solution to the differential equation ( 25 ) . For example, a 



canonical Lagrangian with speed of sound cs — const. — 
1 is just the case s = 0, so that n — 1 and C = 1, and 
Eq. (25) becomes 



XX 



0, 



with general solution 



C = f{^)X-V{<l>). 



(27) 



(28) 



Here / (0) and V (0) are free functions which arise from 
integration of the second-order equation (25). The func- 



tion / ((f)) can be eliminated by a field redefinition dip = 
\/ f (0)d(/>, resulting in a manifestly canonical Lagrangian 
for if, as we would expect from setting cg = 1. We 
emphasize that Eq. (251 is constructed using the so- 
lution (19), and is not a general condition on the La- 



grangian. That is, Eq. ( 25 ) allows us to construct a La- 



grangian which admits solutions of the desired form, but 
those solutions are not necessarily unique. A canonical 
Lagrangian can support inflationary solutions, but not 
tachyacoustic solutions, and is therefore not of interest 
here. However, other choices of n do yield tachyacoustic 
solutions, and we focus on two such choices: 

1. n = 0: A Cuscuton-like model. 

2. n 3: A DBI model. 

In the next section, we consider the attractor behavior of 
the solution given by Eqs. ( |17|18|19| in both the Cuscu- 
ton and DBI cases. 



III. ATTRACTOR BEHAVIOR 



The Cuscuton Case 



The case where the exponent n = in Eq. ( 25 ) corre- 



sponds to taking s = —2s, and the resulting Lagrangian 

is m 



C{X,^)^2f{<j,)^/X + CX-V{4>), (29) 

where / {cf)) is an arbitrary function, V (0) is the potential 
and the constant C in Eq. ( 24 ) is given by 



C - ^^p' 



(30) 



The power-law solution ( 17|18|19" l corresponds to the 
choice 



/(0) 

y(<^) 



V2Ml,Hoe 
s4>o 



'-/s+l 



'0 

2e/s 



,(31) 
(32) 



The Hubble parameter ^ is 



H = 



^02 + ^(0) 



(33) 



and the sound speed ([s]) is 



Cs 




The equation of motion for the Lagrangian ( 29 ) is 

'<j} + iHcl^ + C-^V' (0) =0, 
for which the exact solution ( 17|18|19 ) is 



(34) 



(35) 



= -sH, 



090 



H (0) = i/o 



cs (0) = 



(36) 



It is straightforward to verify that this solution corre- 
sponds to power-law evolution. 



a cx t^^"^ tyi e 



-N 



(37) 



where the Hubble parameter and speed of sound evolve 
according to Eq. (16). The resulting model is tachya- 



coustic, that is, has a decaying, superluminal speed of 
sound and shrinking acoustic horizon in comoving units. 
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It can be shown [7] that such cosmologies produce a 
power-law spectrum of superhorizon curvature perturba- 
tions with spectral index 



1 



2e + s 



1 



(38) 



so that the scale-invariant limit is s = — 2e. A red- 
tilted scalar spectrum index of perturbations such as 
that favored by WMAP, n = 0.96 ± 0.026 ^ can be 
achieved if we choose the value of the flow parameter s 
to lie in the interval s = [-1.907e, -1.980e]. (For the 
remainder of this paper we will consider the scale in- 
variant limit s = — 2e, which is sufficient for the pur- 
pose of demonstrating a dynamical attractor.) In this 
sense, tachyacoustic cosmology represents an interesting 
possible alternative to inflation: instead of accelerating 
expansion, tachyacoustic models produce a spectrum of 
perturbations consistent with the data via a decreasing, 
superluminal sound speed in (for example) a matter- or 
radiation-dominated background. However, in order to 
represent a viable cosmological solution, the analytic so- 



lution (36) must be a stable dynamical attractor, so that 



many different choices of boundary condition will con- 
verge on a single solution at late time. 



To study the attractor properties of the solution ( 36 1 , 
we construct a phase-space representation of the full 
equation of motion ( 35 ) . We introduce dimensionless 
phase space variables 



y 



3s2 i/o0o 



(39) 



The equation of motion can then be written as a phase- 
space evolution equation using |21) 



d(j) 



— ^ — y{x)y (x), 



(40) 



so that 



y (x) y' {x) 3c| 




V [x) = 0. 



(41) 

We define dimensionless versions of the potential V and 
Hubble parameter H as 



V {x) 



V 



.) 



x'^ 



" (l - l^') , (42) 



and 



so that the phase space evolution equation ( 41 1 reduces 
to the simple form 

V {x) y' {x) -f iclh {x) y {x) + v' {x) = 0. (44) 



Tachyacoustic attractor (cusciiton): e = 2, s - 




FIG. 1: Evolution of the dynamical variables x and ?/ in a 
radiation-dominated tachyacoustic model with a cuscuton La- 
grangian. The Hamilton-Jacobi trajectory is the red (dashed) 
line. 



Similarly, we define a dimensionless warp factor 
' 2e s(f>Q 



- V 3 721/2 i/oe 
The sound speed is then 



4 



1 



y{x)' 



(l-x2). (45) 



(46) 



In terms of the dimensionless variables, the analytic so- 
lution ( 36 1 becomes 




/s+l 



(47) 



We evaluate the equation of motion (44 1 numerically 



to demonstrate that the solution (47) is a phase-space 
attractor. Figure [T] shows attractor behavior in the 
radiation-dominated case (e = 2), and Figure [2] shows at- 
tractor behavior in the matter-dominated case (e — 3/2). 



B. The DBI Case 



Similar to the cuscuton case considered in the previous 
section, we can construct a Dirac-Born-Infeld Lagrangian 
by choosing s = s, so that n = 3 in Eq. (25), which 
corresponds to the Lagrangian [7] 



c {X, 0) = -r 1 (0) ^i-f{(t>)x + r 1 (0) - y (0) . 

(48) 

For e and s constant we obtain an exactly solvable sys- 
tem, with 



-4e/s 



1 - 



1 



(49) 
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Tachyacoustic attractor (cuscuton): e = 1.5, s = -3. 



FIG. 2: Evolution of the dynamical variables x and ?/ in a 
matter-dominated tachyacoustic model with a cuscuton La- 
grangian. The Hamilton- Jacobi trajectory is the red (dashed) 
line. 



and 



/ 1 



4c/s-2 



The sound speed ^ for the DBI Lagrangian is 



cs 



c 



X 



(50) 



(51) 



and the Hubble parameter ^ is 
1 



3A/2 
1 

3A/| 



i2XCx - C) 



l~cs 



The equation of motion for for the field (j) is then 



with solution 



s 

2 \(i)o 



-2£/s+l 



r(0) 



(52) 



= 0, 
(53) 



H {<!>) = Ho 



Cs (0) 



(54) 



This solution again corresponds to power-law evolution 
a oc oc t^/*^. As in the cuscuton case, the spectral 
index of scalar perturbations is 



n = 1 



2e + s 



(55) 



l-e-s' 

so that the scale-invariant limit corresponds to s = — 2e. 



We study the attractor properties of this solution by 
defining dimensionless variables 



(56) 



/3HnM, 



and 



1 

1 + a;2 



g{x) EE iMlHlf (0) = (1 - x') (57) 



The dimensionless Hubble parameter is 



h(x) ^ - , 

VSMpHo V 

where the sound speed is 



8e 
3s2 



1 - Cs (x) 
cs (x) g {x) 



+ V {x) 



(58) 



(59) 



In terms of these dimensionless variables, DBI equation 
of motion is: 



y (x) y' ix)+3hcly ix) + ^^-^y' (^) + (4-l)4+4^;' = 0. 



25 



(60) 



The analytic solution (|54|) is 



h{x) = 



Cs ^x 



2 V_3s2 
Be 




2e/s 



(61) 



We evaluate the attractor behavior of this solution in 
the scale-invariant limit s — — 2e. Figure ^ shows the 
radiation-dominated case (e = 2), and Figure Q shows 
the matter-dominated case (e = 3/2). In both cases, the 
analytic solution (54 1 is a dynamical attractor. 



IV. CONCLUSIONS 

In this paper, we consider the dynamical stability of 
"tachyacoustic" cosmological models [7], which generate 
superhorizon cosmological perturbations via a decreas- 
ing, superluminal sound speed instead of accelerating ex- 
pansion (as in the case of inflation). It is known that such 
cosmologies can produce nearly scale-invariant scalar per- 
turbations, consistent with current data. However, it has 
not been previously shown that such solutions also cor- 
respond to dynamical attractors, which is a necessary 
condition for such models to be cosmologically viable. 
Such models are realized in scalar field theory via a non- 
canonical Lagrangian [S]. In our analysis, we have con- 
sidered two particular choices of Lagrangian which give 
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Tachyacoustic attractor (DBI): e = 2., s = -4. 




X 



FIG. 3: Evolution of the dynamical variables x and y in 
a radiation-dominated tachyacoustic model with a DBI La- 
grangian. The Hamilton- Jacobi trajectory is the red (dashed) 
line. 



Tachyacoustic attractor (DBI): e= 1.5. s = -3. 




speed, resulting in a scale-invariant primordial power 
spectrum. The first case is a so-called "cuscuton" La- 
grangian, which is linear in the field kinetic term instead 
of quadratic as in the case of a canonical Lagrangian. 
We numerically integrate the full phase space for the field 
evolution, and show that the power-law tachyacoustic so- 
lution is in fact a dynamical attractor. Such models are 
of particular interest because they predict a detectable 
contribution to cosmological non-Gaussianity |15j . Sec- 
ond, we consider a Dirac-Born-Infeld (DBI) Lagrangian 
giving identical power-law behavior and show that the 
power-law solution to this Lagrangian is also a dynami- 
cal attractor in the full phase space. We present results 
for the scale-invariant limit, but considering a slightly 
"red" spectrum as favored by data does not alter the at- 
tractor properties of the solution. We conclude that, like 
inflation, tachyacoustic cosmology can generate a scale- 
invariant power spectrum via dynamically stable cosmo- 
logical evolution. We note that in this analysis we con- 
sider only the classical stability of tachyacoustic models. 
This may not apply at the quantum level. Negative- 
tension branes, for example, are known to have instabil- 
ities at the quantum level [HI [53]. However, we are not 
aware of any general argument indicating that theories 
with a superluminal sound speed are necessarily unstable 
to quantum mechanical flucutations. 
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